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We derive generic relativistic hydrodynamical equations with dissipative effects from the un- 
derlying Boltzmann equation in a mechanical and systematic way on the basis of so called the 
renormalization-group method. A macroscopic frame vector is introduced to specify the frame on 
which the macroscopic dynamics is described. Our method is so mechanical with only few ansatz 
that our method give a microscopic foundation of the available hydrodynamical equations, and 
also can be applied to make a reduction of the kinetic equations other than the simple Boltzmann 
' equation. 
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The relativistic hydrodynamical equation is widely used in various fields of physics, especially in high-energy nuclear 
physics and astrophysics [l], 0, E| , and it seems that the study of the dissipative hydrodynamical equations is now 
■ becoming a central interest in these fields 0, H i, @, 0, 0, B G3 • 

It is, however, noteworthy that we have not necessarily reached a full understanding of the theory of relativistic 
hydrodynamics for viscous fluids: (1) It is well known that the relativistic hydrodynamical equation with dissipation 
has a form depending on the choice of the Lorentz frame or the definition of the hydrodynamical flow. The typical 
• frame is that on the particle or the energy flow, and the typical corresponding equations are the ones of Eckart 
PH and Landau [l2j], which were derived phenomenologically on the basis of the number and energy-momentum 
conservation laws, the second law of thermodynamics and some specific assumptions on the choice of the flow. One 
might consider that the difference of the equations should be merely a kind of choice of coordinate to describe dynamics 
in a easy and practical way, and one might be able to go back and forth between the two reference frames, by a Lorentz 
transformation. However, it may not be the case. To be more specific, let ST' 11 ' and SN^, be the dissipative part of the 
*^ energy-momentum tensor and the particle current, respectively. The point is that they are not determined uniquely 
only by the second law of thermodynamics without so me p hysical constraints involving the four flow-vector it M with 
Q_ij u^u^ = 1: The constraints of Eckart are expressed as [LI, (i) u li u„ST ,iV = 0, (ii) u^SN^ — 0, (iii) A^SN" = 0, 
while those of Landau are (i) u /x u„5T' iI/ = 0, (ii) u^SN^ = and (iv) u f _ l A up ST fJ - 1 ' = 0, where A pv = — u^u^. 
As one sees, the first two constraints are the same for the two theories, and the respective third condition specifies 
J> , the frame of the coordinate. It is here noteworthy that there is a proposal by Stewart for the constraints for the 
Eckart frame, i.e., the particle-flow frame, as given by (v) ST^ = 0, (ii) u^&N^ = 0, (iii) A^SN 1 ' = 0, where the 
condition (i) of Eckart is replaced by a different one (v). One may ask if both the Eckart and Stewart constraints 
make sense or not. It has been scarcely examined [l5[ whether the constraints proposed phenomenologically hold 
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internal consistency or match the underlying microscopic kinetic theories. (2) It is also known that these equations 
containing the time derivative only in the first order, unfortunately, suffer from the unphysical instability [H, fl6| ; 
the instability is attributed to the lack of the causality, and invented is a fluid dynamical equation which contains 
the second-order time-derivative as well as the first-order one [ill, E3, G3] • The new equation called Israel-Stewart 
equation [2(| contains new transport coefficients like the relaxation time yet to be determined microscopically or 
phenomenologically 0, GJ, EI E0| • 

There are some attempts to derive the phenomenological equations from the microscopic equation, i.e., the relativis- 
tic Boltzmann equation, with use of the Chapman-Enskog method [2l| and the Grad's fourteen- moment method [22| . 
and hence give a microscopic foundation to them 0, [IH, QjJ He], [HI ■ We notice that such a microscopic approach has 
another merit that the hydrodynamical equations can be obtained, which are consistent with the underlying kinetic 
equation, because the theory gives the theory connect the kinetic and the hydrodynamical regimes and hence gives 
the natural initial condition for the hydrodynamical equation. This point should have a practical merit when one 
analyzes the freeze-out regime and also the hadron corona in the expanding system such as the intermediate stage of 
the relativistic heavy- ion collision 0, H, H3| • 

One must say, however, that the past works in the microscopic approaches are not fully satisfactory in giving the 
foundation to the phenomenological equations: a. Although the past works succeeded in identifying the assumptions 
and/or approximations to reproduce the known hydrodynamical equations by Eckart, Landau and Israel, one must 
say that the physical meaning of these assumptions/approximations is somewhat obscure, and so is the uniqueness 
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of those hydrodynamical equations as the long-wavelength and low-frequency limit of the microscopic theory, b. It is 
also to be noted that the Chapman-Enskog method and the Grad's fourteen-moment method themselves have some 
ad-hoc part and may not be fully mechanical and systematic methods for the reduction of dynamics. 

The purpose of this article is to derive the relativistic hydrodynamical equations with dissipative effects from 
the microscopic theory in a more natural and systematic way, and thereby establish the microscopic foundation of 
them. We shall show that when the particle-flow frame (Eckart frame) is chosen, the constraints by Stewart but 
not Eckart must be taken, while the resulting hydrodynamical equation which manifestly satisfies the second law of 
thermodynamics is different form that given by Stewart. 

The problem is how to reduce the dynamics to a slower dynamics using a reliable reduction theory of the dynamics 
[25l ] . which should be as mechanical as possible. In fact, van Kampen [23| applied his method [26[ of the reduction 
to derive the relativistic hydrodynamic equations for the viscous fluid. However, his method is not formulated in a 
covariant way and hence failed in reproducing the phenomenological hydrodynamical equations with dissipation. We 
shall a ppl y another powerful reduction theory known as the "renormalization-group(RG) method" [13, HI, [2f| l30l . 
I3TL [32l. l33| to reduce the relativistic Boltzmann equation to the hydrodynamical equations. Our theory is manifestly 
covariant and relies on almost no assumptions. 

The RG method is a systematic reduction theory of the dynamics leading to the coarse-graining of temporal and 
spatial scales. This method is suitable for our present purpose to elucidate the physical meanings of each process 
of the reduction since the method is so mechanical and does not require any specific assumptions. Indeed, the RG 
method is already applied satisfactorily to the derivation of the Navier-Stokes equation from the (non-relativistic) 
Boltzmann equation [32I [331 ] . The present work is an attempt to extend these previous works to the relativistic case. 

In this article, we shall treat, as a simplest example, the classical but relativistic system composed of identical 
particles; an extension to multi-component systems is possible and will be discussed in a separate paper [34[. Then 
the relativistic Boltzmann equation [18| describing such a system reads 




where f P (x) denotes the one-particle distribution function defined in the phase space (x M , p M ) with p M being the 
four momentum of the on-shell particle; p^p^ = p 2 — m 2 and p° > 0. The right-hand side of Eq.{T]) is called the 
collision integral, C[f] p (x) = ^ J2 P2 ^ J2 P3 Pi\P2 , Pa) (f P2 (x) f P3 (x) - f p (x) f Pl (x)) , where uj(p,px\p 2 , p 3 ) 

denotes the transition probability owing to the microscopic two-particle interaction and contains the delta functions 
5 4 (p + pi — pi — Pz) representing the energy-momentum conservation. The transition probability has the symmetric 
properties due to the indistinguishability of the particles and the time reversal invariance in the scattering process: 
lj(j>, pi\p2 , Pa) = u(p 2 , P3\p, Pi) = w(pi , p\p3 , pa) = Lu(p 3 , p 2 |pi , p). It should be stressed here that we have 
confined ourselves to the case in which the particle number is conserved in the collision process. 

The particle-number and energy-momentum conservation in the collision process ensures the following important 
relations; J2 P ^ C[/]p( x ) = S P jp P 11 C[/]p( x ) — 0, which says that the five quantities (1 , p^) are collision invariants. 
We have the continuity or balance equations for the particle current and the energy-momentum tensor, 



respectively. Notice that while these equations have the same forms as the hydrodynamical equations, nothing about 
the dynamical properties is contained in these equations before the evolution of the distribution function f p (x) is 
obtained by solving the microscopic Eq.|T]). 

The entropy current may be defined by = — J2 P pTS P^ fp( x ) (^ n fp( x ) ~ 1); which is conserved only if lnf p (x) 
is a collision invariant, or a linear combination of the basic collision invariants (I , p M ) as hif p (x) = a(x) + p^ ^(x) 
with a(x) and /3^(x) being arbitrary functions. 

Since we are interested in the hydrodynamical regime where the time- and space-dependence of the physical quan- 
tities are small, we solve Eq.([T]) under the situation where the time variation of f P (x) is slow with long wavelengths. 

To retain the Lorentz covariance, and achieve a coarse-graining, it is found convenient to introduce a macroscopic 
Lorentz vector ap 1 characterizing the flow, which specifies the covariant coordinate system (frame) and we call the 
macroscopic-frame vector. Although af^ may depend on the momentum p and the space-time coordinate x, the time 
variation of it is supposed to be much smaller than that of the microscopic processes. Our point is that the separation 
of the scales can be nicely achieved by the RG method. 

Keeping in mind the above scale difference, we define the covariant coordinate system (r , <r M ) as dr = a^(x) dx^ and 
e' 1 da^ = [g^ — a^(x) a p {x) / 'ai{x)) dx v = A^(x) dx v . We notice that the small quantity e have been introduced 
to tag that the space derivatives are small for the system we are interested in; the e may be identified with the ratio 




p"0 M / p (s) = C[f] p (x), 
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of the average particle distance over the mean free path. In this coordinate system, Eq.Q]) reads 
d 1 1 

fl-M T . a ) = 7 r^/p T, cr) -E rP'V/ p T, cr), (3) 

where a^(T, cr) = a^(x), A^(r, cr) = A^ u (x) and / p (r, cr) = f v (x). Since £ appears in front of V M = 
A^ ly (r, cr) g^77, Eq.Q has a form to which the perturbative expansion can be possible. In fact, a natural but 
significant assumption underlies this seemingly mere rewrite of the equation that only the spatial inhomogeneity over 
distances of the order of the mean free path is the origin of the dissipation. We notice that the RG method applied 
to non-relativistic Boltzmann equation with the corresponding assumption nicely leads to the Navier-Stokes equation 
[33. l33l] : the present approach is simply a covariantization of the non-relativistic case. 

Some remarks on a^(r, cr) are in order: (i) In the perturbative expansion, we shall take the coordinate system 
where a^(r, a) has no r dependence, i.e., a^(r , cr) = a^(cr). (ii) In our coordinate system, the physical quantity 
transported per unit r is given by p ■ a p (a), we can control what is the flow represented in our theory by varying the 
specific expression of a^(cr). Owing to this freedom inherent in our coordinate system our theory may lead to various 
hydrodynamical equations, including the ones in the Eckart and Landau frames for non-ideal fluids. 

In accordance with the general formulation of the RG method given in [2^, [29|, [3l|, we first try to obtain the 
perturbative solution f p around the arbitrary initial time r = tq with the initial value / p (to , cr); f p (r = tq , cr ; tq) = 
fp(j~o , cr), where we have made explicit that the solution have the To dependence. The initial value is not yet specified, 
we suppose that the initial value is on an exact solution with To being varied. The initial value as well as the solution 

are expanded with respect to e as follows; f p (r , a; t ) = /p 0) (r , a ; r )+e/( 1) (r,cr;r )+ e 2 /i 2) (T , a; t )-\ , and 

fp( T o 7 cr) = fp " 1 ( T o 1 cr) + e fp 1 ^ (tq , cr) + £ 2 fj^ (tq , cr) + • • • . The respective initial conditions at r = To are set up as 
fp^ (to , cr; To) = fp (tq , a) for I = 0, 1, 2, • • • . In the expansion, the zeroth-order value f p °^ (tq , a ; To) = f p °^ (to , cr) 
is supposed to be as close as possible to an exact solution. 

Substituting the above expansions into Eq.Q in the T-independent but To-dependent coordinate system with 
a p( T j a ) ~^ a pi T o , f) = a p (cr; t ), we obtain the series of the perturbative equations with respect to e. The 
zeroth-order equation reads 

Since we are interested in the slow motion which may be realized asymptotically when r — > oo, we should take 
the following stationary solution or the fixed point, ^r/p°' ) = 0, which is realized when C[f^] p = for arbitrary 

cr. This implies that ln/p ' is a linear combination of the five collision invariants (1 , j/ 1 ), and hence /p ' is found 
to be a local equilibrium distribution function or the Juettner function (the relativistic analog of the Maxwellian) : 
fp°\T , cr ; t ) = (2tt)" 3 exp [(/i(cr ; t ) - p^ u p (a ; to)) / T(a ; t )] = f p q (a ; r ) with u^(cj ; r ) u^(a ; t ) = 1, which 

implies that fp ^ (r , cr ; t ) = (tq , cr ; t ) = (t , cr) = /p q (cr ; To). It should be noticed that the five would-be 
integration constants T(a; To), ^(cr; To) and u M (cr ; to) is independent of t but may depend on To as well as a. 
Next the first-order equation reads 

where the linear evolution operator A and the inhomogeneous term i* 1 are defined as A pq = (p ■ a p ) _1 ^-C[/] p |^ _^ cq 

and F p = —(p- a p ) _1 p • V/p q , respectively. To obtain the solution which describes the slow motion, it is convenient 
to first analyze the properties of A, especially its spectra. For this purpose, we convert A to another linear evolution 
operator L = / oq_1 Af cq with the diagonal matrix f p q = fp q S pq . Let us define the inner product between arbitrary 
vectors ip and ip by 

( V ,^)=Y, \(p-ap)P q V p^ p . (6) 
p P 

With this inner product, one can nicely show that L becomes self-adjoint {if, Lip) = (L(p,ip) and non-positive 
definite (ip , L(p) < 0, which means that the eigen values of L are zero or negative. The eigen vectors of the zero 
eigen value are ipQ p = p^ for a = [i = ~ 3 and <p\ p = m, which span the kernel of L and satisfy L (p^ — 0. Notice 
that these five vectors are the collision invariants. Following [3l|, we define the projection operator P onto the kernel 
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of L which is called the P space and the projection operator Q onto the Q space complement to the P space: 

[PV] p = $>^<^>, (7) 

a/3 

with Q = 1 — P where rf^a is the inverse matrix of rj a ^ = (<Pq , <Pq)- 

In the following, we shall suppress the argument a and the subscript p when no misunderstanding is expected. The 
solution to Eq.© with the initial condition / (1) (r = r ; r ) = / (1) ( T o) is found to be / (1) (V ; r ) = e (r " ro)A {/ (1) (r ) + 
A' 1 QF} + (t-t )PF- A' 1 Q F, where the modified projection operators P = f cq P / cq_1 and Q = f eti Q / cq_1 have 
been introduced. Notice that the first term would be a fast motion coming from the Q space, which we can nicely 
eliminate by choosing the initial value /^( r o)j which has not yet been specified; thus we have the first-order solution 

fW(T;T ) = (T-T )PF-A- 1 QF (8) 

with the initial value determined as /^(t ; r ) = /^( T o) = —A' 1 Q F. We notice the appearance of the secular 
term proportional to r — To, which apparently invalidates the perturbative solution when |t — to| becomes large. 
The second-order equation [3|| is written as 

l_p) =A ~ f m + { r- TQ )H + I, (9) 

with H p = —(p ■ dp)' 1 p ■ V[PF] P and I p = (p ■ a p ) _1 p • V[A _1 Q F] p , the solution to which is found to be 

/(2)(t;t ) = "^{f^i^+A-'QH + A-'Ql} 

+ \ (r -t ) 2 PF+(t-t ){PI- A' 1 QH}- {A~ 2 QH + A' 1 Ql}. (10) 

Again the would-be fast motion can be eliminated by the choice of the initial value /^(tq), and thus we have the 
second-order solution 

f {2) (r ■ t ) = i (r - T ) 2 PF + (r - r ) {P I - A' 1 Q H) - {A~ 2 QH + A^Ql}, (11) 

with the initial value determined as /^(to ; to) = f^(ro) = —A~ 2 Q H — A' 1 Ql. We notice again the appearance 
of the secular terms. 

Summing up the perturbative solutions up to the second order, we have an approximate solution around t ~ To to 
this order; /(t ; t ) = /^°'( T ! T o)+£ f^\ T ', T o)+e 2 f^Hr ; T ) + O(e 3 ), containing the secular terms which apparently 
invalidates the perturbative expansion for t away from the initial time t . 

The point of the RG method lies in the fact that we can utilize the secular terms to obtain a solution valid in a 
global domain. Now we may see that we have a family of curves /(t ; To) parameterized with To. They are all on the 
exact solution f(r) at r = tq up to 0(e 3 ), but only valid locally for t near To- So it is conceivable that the envelope E 
of the family of curves which contacts with each local solution at r = To will give a global solution in our asymptotic 
situation. According to the classical theory of envelopes, the envelope which contact with any curve in the family at 
t = t is obtained by 4rJv( T > a 5 T o) | To=T = °j or explicitly 

£{f? - e [A" 1 QF] p } - e [PF] p s 2 { [Pl] p [A+QH],} + 0( £ 3 ) = 0. (12) 

This envelope equation is the basic equation in the RG method and gives the equation of motion governing the 
dynamics of the five slow variables T(t), h{t) and it^(r) in /° q . The global solution in the asymptotic region is given 
as the envelope function, 

f p (T, a ■ to = t) = f p (r , a) = / p cq -s[A- 1 QF] p - e 2 {[A~ 2 QH] p + [A' 1 Ql] p } + 0(e 3 ), (13) 

where the exact solution of Eq. (ll2|) is inserted. Thus one sees that / p (t) now describes the macroscopic-time evolution 
of the one-particle distribution function in Eq.([T]), because the time-derivatives of the quantities in f p (r) are all in 
the order of e. We emphasize that we have derived the infrared effective theory of Eq.{T]) in the form of the pair of 
Eq.'s (fl"2"l) and (TIB")) . This is one of the main results in this article. 
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Now the RG /Envelop equation Eq. (fT2| is actually the hydrodynamic equation governing the the five slow variables 
T(t), and m m (t). To show this explicitly, we apply P from the left and then take the inner product with the five 
zero modes ipfi [37|. Putting back e = 1, we arrive at d^J^ dro = with 

■Cro = E/"*{/; q -[ i " 1 o f L} = ^+^ ( 14 ) 

p 

where Jg" = J2 P jp P^ <Po P fp q an d SJ^ a = — J2 p ^P^fop [A~ 1 QF] p - Here Jq" represents the currents of perfect- 
fluid part, while SJ^ a the dissipative part. In deriving Eq. (fT4|) , we have used the following relation : ~^2 p (p ■ 
a pX P [ p ^] = Ep jo (P • Op) Pop ^P) obtained from the definition P = /eqp/eq-i and Eq.©. Notice that the 
derivation of Eq. (|14[) is accomplished for arbitrary a p which is now dependent on r as well as on a; in other words, 
we have a set of relativistic hydrodynamic equations for viscous fluids, which have still the freedom of choice of ap 1 . 
The is also one of the main results in the present work. 

In the rest of the paper, we shall show how known relativistic dissipative hydrodynamic equations are derived with 
a choice of the macroscopic frame vector ap 1 . 

As a simple but nontrivial choice, let us take the following set of the macroscopic frame vectors with 9 being a 
constant; = (p-u)' 1 {(p-u) cos 0+m sin 9} =6%, which leads to the following relations A. p iV — g^ v —uP 1 u v = A^ v ; 
then one finds that (p • u)' 1 {(p • u) cos9 + m sin6>} ^ = dp = D and V p = d v = V^ 1 . We notice that the 
factor m is introduced simply to make the expression dimensionless, so our method is also applicable to the case of 
massless particles. 

In the reduction of the currents in terms of the fluid dynamical quantities, the central task is the evaluation of 
A' 1 Q F. A straightforward but tedious manipulation leads to the following formula, 

[A- 1 QF] p = f? J2 °v\ \ ( n « X <> + % X ^ + nr *v) . (15) 

where C pq = {p ■ 9 p )L pq being independent of 8^. Here we have introduced the following quantities: the dissipative 
currents, i.e., the bulk pressure, = (4/3 — 7) (p • u) 2 + [(7 — 1) T h — 7 T] (p ■ u) — 1/3 m 2 , the heat flow, = 
-\(p-u)-Th] A» v p v and the stress tensor, H^ = 1/2 (A pp A 1 " 7 + A P<T A vp - 2/3 A pv A P<T ) pp p a . The coefficient of 
each dissipative current is the respective dissipative thermodynamic force, i.e., Xg = (—z 2 cos 2 0) V • u/{z 2 cos 2 9 + 
z 2 (3 7 - 4) sin 2 (9 -3z[l-(h- 1) (7 - 1)] cosfl sin<9}, X 9 p = V^lnT - cos 9 V M ln(nT)/(h cos(9 + z sin0), and 
Xp U = 1/2 (Ap P A va + Ap a A vp — 2/3 Ap V A pa ) V p u CT , where the following definitions have been used; the reduced 
mass z = Tjr, the particle density n = (2 n)~ 3 47rm 3 eT {z _1 Kz(z)}, the energy density e = mn{Ks(z) / K2(z) — z -1 }, 

the reduced enthalpy per particle h = e +y T and the ratio of the heat capacities 7 = 1 + {z 2 + 3h — (h— l) 2 }" 1 , with 
the modified Bessel functions K 2 (z) and K 3 (z). 

Since the dissipative part of the currents has been obtained for arbitrary 9, we can obtain the microscopic formulae 
for the transport coefficients, i.e., the bulk viscosity C, the heat conductivity A and the shear viscosity 77, which are 
nicely given in forms of the Kubo formula as follows: Using the new inner product ( ip , i/j ) eq = ^ p ^ f p q (p p tp p and the 

evolved microscopic dissipative currents [n(s)] = J2 q [ eS C ] pq anc ^ so on ' wc navc C = {^-/T) J ° O ds( n(0) , U(s) ) cq , 
A = -(1/3T 2 ) ^^(^^(O) , Q^s) ) and 77 = (1/10T) J °°ds( n p,y (0) , ^^(s) ) It is noteworthy that these trans- 
port coefficients are independent of 9 while the dissipative thermodynamic forces are dependent on it. 

Although we can now write down the generic form of the currents jQ ro for arbitrary 9, we shall give here them 
only for a few values of 9 for the sake of space, but show that our generic currents included the known currents for 
the relativistic fluid dynamics for a viscous fluid. 

(A) With the simplest choice a p = i.e., 9 = 0, we have 

( eu^u u - (p + (X)A^ + 2rjX^ a-v, 

J ^° = \mnu»-Z\X» a = 4, (16) 

I h 

where p = nT, X = —V • u and X^ 1 = V M lnT — h' 1 V M ln(nT). If we employ the constitutive equations; H = ( X, 
Q p = T XX P and Ii pv = 2rjX^ L ' — HA* 1 ", it is possible to rewrite the equation in terms of the flows rather than the 
thermodynamic forces as 

( eu^u v -pA^ + W U a = v, 

J ^ d ™ = I mnu^—^—Q^ a = 4. 

I hT 
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One should notice that J^ydro m Eq. fTBf or Eq. (fT7|) agrees completely with the Landau theory for non- ideal fluids. 
This was actually anticipated: In fact, if we take a natural choice that = it^, the physical quantity transported 

per unitr becomes the microscopic thermal energy (p ■ u), which is reduced to the familiar form ii| 2 in 
the non-relativistic limit. Thus it is natural that the currents JyLdro w itfi the choice ap 1 = becomes those in the 
Landau frame. 

(B) Another simple choice for the macroscopic frame vector is given with 9 = tt/2; a£ = j^u^. A similar 
manipulation as the above gives the following currents 

f U + ZC,X)u^u v - (p+CX)A tlv + XTu^ X v + XTu v X^ +2nX tlv a = v. 

•Cdro = \ * , ( 18 ) 

I mnu^ a = 4, 

where X = -{1/3 (4/3 - 7)" 1 } 2 V • u and X^ = V^ 1 InT. In terms of the flows, we have 

f eu»u v -pA» v + ZIlu»u v + u tt Q v + u v Q 11 + a = v, 

J Cdro= , , (19) 

y mnu^ a = 4, 

where II = ( X , — XT X^ and IF" = 2r)X^ v -ft A"". Since the dissipative term with the heat conductivity is 
absent in J^t dm , one might be tempted to identify the above equation with the Eckart one, but it is not the case: A 



detailed examination tells us that Eq. (fTg|) is rather a modified or corrected version, we would say, of Stewart's[14|. 
This observation is based on the fact that Eq. (fl~8|) meets the Stewart's constraints, (v), (ii) and (iii) but not the 
Eckart's, (i), (ii) and (iii), mentioned in the Introduction. The reason why we call Eq. (fl~8|) the corrected version 
of Stewart equation will be shortly explained. A noteworthy point is the following: In the Eckart theory [38] used 
conventionally, the appearance of the terms with the bulk viscosity is the same as the Landau theory. In contrast, 
Eq. (|18p shows the different appearance so that the bulk viscosity has an influence on the terms of not only the pressure 
p but also the energy density e. We can show that our equation (fT5|) satisfies the second law of thermodynamics as 
will be explicitly shown in a separate paper [34L l35j|. although the original Stewart equation may not. 
(C) The other interesting choice of 9 is 6 = — 7r/4, for which we have a>i — (P'")~ m u m ana ; 

y ' ' P \J2p-u 

^ j tu^u v -pA^ - XTu^ X v ~ XTu v X^ + 2f 1 X> nJ a = v, 

^hydro = ) u . _ -ry u . (^0) 

[mn/-ATF a = 4, 

where X^ = z(h — z) _1 V M lnT — z(h — z)~ 2 V M ln(n T). The alternative expression in terms of the flows is 

eu^u" -pA» u -u^Q" -tt^Q^ + n^ a = v, 
mnu^ — Q M a 



^hydro — { „ ^„ , ( 21 ) 



where = XT X^ and IP 1 ' = 2r\X liV . In this frame, X = Xg = _ n u — so that the associated bulk pressure 
ff(= (X) disappears. The unique feature of this equation is the absence of terms with the bulk viscosity. As far as 
we know, this type of the dissipative fluid dynamical equation is written down for the first time, which was made 
possible by the introduction of the macroscopic frame vector in the powerful RG method. 

Now we have seen that the fluid dynamical equation for dissipative fluids in the Eckart frame does not take the 
form of Eckart, if we start from the underlying Boltzmann equation. This fact is actually known for some time [l5l |. 
Nevertheless it is instructive to see that there can not exist the frame vector that leads to the Eckart's constraints, 
(i), (ii) and (hi), because it would help to elucidate the physical meaning of the P- and Q-space in our formalism and 
their relation to the phenomenological constraints imposed to the dissipative part of the energy- momentum tensor SJ^" 
and the particle-number current SJ^ 4 , which are defined just below Eq. (fT4"|) . Eq. (fl~5|) tells us that the dissipative part 
of the one-particle distribution function f p , i.e., — [A -1 Q .F] p can be written as fp q <pp with (f> p = — [L' 1 Q / cq_1 F~\ . 
Since the Q-space vector (f> p is orthogonal to the five P-space vectors ipQ p by definition, we have 

(pft, <j>) = for a = 0, 1, 2, 3, 4. (22) 

Here the inner product is defined by Eq.((6]) where ap 1 enters with the form (p ■ a p ). The notable point is that these 
five identical equations exactly correspond to the constraints on 8J^ V and 5J^ A . In fact, Eq. (|2"2"|) with the choice 
of (A) — gives J2 P (p ' u ) fp q '■Pop §v = 0, which is equivalent to the set of equations u u 8J IJ ' v — and 
SJ^ = 0, because 5 J^ a — J2 P P^ fo p fp 4>p- The former equation may be re-expressed by u v 8J^ V = and 
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A MP u„ (5J A " / = 0. Thus one can readily see that these equations coincide with the Landau's constraints, (i), (ii) and 
(iv). Similarly, with the choice of (B) — ^ u M , we have m ifQ p <fi p — 0, which means that S J^ 4 = and 

SJ^ = 0; the former equation is equivalent to the set of equations 5 J m4 = and A M „ 5 J^ 4 = 0. Thus one sees 
that the Stewart's constraints, (v), (ii) and (iii) are derived. Here we have used the on-shell condition m 2 = p^p^,. 
Now it is easy to show that there exists no leading to the Eckart's constraints, (i), (ii) and (iii), simultaneously. 
In order to lead to the constraints (ii) and (iii) on SJ^ 4 , Eq. lj2"2"|) with a = /i requires that (p ■ a p ) = const., i.e., 
independent of p^. On the other hand, in order to lead to the constraint (i) on SJ^ V , Eq. ([22|) with a = 4 must lead 
to (p ■ a p ) — const, x (p • w) 2 , which is in contradiction with the condition for (ii) and (iii). Thus we can conclude 
that there exists no ap 1 leading to the Eckart's constraints (i), (ii) and (iii) simultaneously. 

What we have seen is that whenever the Eckart frame is taken where the particle flow is constructed so as to have 
no dissipative part, the constraint on the dissipative part of the energy- momentum tensor must satisfy (v) but not (i). 
This fact suggests that the Eckart's theory may not be realized as the slow dynamics of the underlying Boltzmann 
equation, and hence nor have microscopic foundation [l5j . 

One of the unique and important feature of the hydrodynamical equations obtained in the present work lies in 
the fact that the thermal force driving the heat flow contains V M T but no Du 11 as in the Landau equation, which in 
contrast with the Stewart and Eckart equations. We shall briefly argue that the Du^ 1 terms should not exist on the 
general ground on the basis of a natural assumption on the origin of the dissipation; the detailed argument will be 
presented in separate papers [3o| . 

First of all, we can trace back the absence of Dw M terms in the thermal forces to our ansatz that the dissipation comes 
solely from the spatial inhomogeneity at the rest frame. This ansatz is clearly expressed in the choice of ap 1 : Although 
the generic form of which keeps the Lorentz covariance can be given as ap 1 = a(p ■ u)u^ + f3(p ■ u) p Vl where 
a{x) and (3{x) are the arbitrary Lorentz-scalar functions, wc have intentionally put f3(x) = so that the differential 
operator V M in the last term in Eq.Q, which gives rise to the deviation from the local equilibrium and hence the 
dissipation, agrees with the spatial derivative, i.e., V M = V M . As mentioned before, the non-relativistic Navier-Stokes 
equation can be nicely derived by the same ansatz that the spatial inhomogeneity leads to the dissipation (32l . |33| , 
and the present approach is simply a covariantization of the non-relativistic case. On the other hand, the ansatz used 
by Landau or Stewart is that the dissipation comes from both of the spatial and the temporal inhomogeneity at the 
rest frame, although no Du^ terms appear in the Landau equation. 

In a separate paper, we shall show that the covariant dissipative hydrodynamic equations constructed in the 
phcnomenological way as was done by Eckart, Landau, Stewart so on but based on the same ansatz as that in the 
present work, completely agree with the equations derived in the present RG method for the Landau and Eckart 
frame, i.e., with the two choices for a£J. 

In summary, we have derived generic covariant hydrodynamical equations for viscous fluids as a reduction of 
dynamics from relativistic Boltzmann equation in a mechanical way with no heuristic arguments on the basis of so 
called the renormalization-group method: This was made possible by introducing the macroscopic frame vector which 
defines the macroscopic Lorenz frame in which the slow dynamics is described. Our generic equation includes not 
only the Landau equationpjj but also a modified Eckart/Stewart equationpjl, [3| and a novel one which does not 
contain the bulk viscosity term. 

We notice that the derived equation is consistent with the underlying kinetic equation, so the equations and also 
the method developed here may be useful for the analysis of the system where the proper dynamics describing the 
system changes from the hydrodynamic to kinetic regime, as in the system near the freeze-out region in the RHIC 
experiment. 

Our method is so mechanical and simple that it is successfully applicable to derive the dissipative hydrodynamics 
for the multi-component system in Landau frame [34| . In the present work, the fluid dynamical variables correspond 
to the conserved quantities in the collision process and span the five-dimensional P space, which is an invariant 
manifold in the terminology of the dynamical systems. One may suspect that the method could apply to derive the 
Israel- Stewart equation [201 ] with fourteen variables in accordance with Grad's moment method [22j . The answer is 
affirmative: We have found (34j that a simple extension of the P-space in our formalism gives the equation equivalent 
to Grad's fourteen- moment method [22| . These extensions will be reported elsewhere. In passing, we emphasize that 
our method itself has a universal nature and can be applied to derive a slow dynamics from other kinetic equations 
than the simple Boltzmann equation. 
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